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Abstract : We construct positive solutions of the semilinear elliptic problem Au + 
Xu + uP = with Dirichet boundary conditions, in a bounded smooth domain O C 
(A^ > 4), when the exponent p is supercritical and close enough to and the parameter 
A G M is small enough. As p ^ solutions have multiple blow up at finitely 

many points which are the critical points of a function whose definition involves Green's 
function. Our result extends the result of Del Pino, Dolbeault and Musso [5] when O is a 
ball and the solutions are radially symmetric. 
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1 Introduction 

In this paper we consider the semilinear elliptic problem 

( Au + Xu + uP = in O 

u > in n (1) 
u = on do, 

where is a bounded regular domain in R^, A > 4, the parameter A S M and the 
exponent p is larger than 

A + 2 

the critical Sobolev exponent. 

When p = Pni Brezis and Nirenberg [3] have proved that (1) admits a solution provided 
< A is less than the first eigenvalue of the Laplacian on with Dirichlet boundary 
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condition. Direct application of Pohozaev's identity [12] shows that solutions of (1) do not 
exist when X < 0, p> pn and 47 is a star-shaped domain. 

In this paper, we are interested in the existence of solutions of (1) in the case where p 
is larger than the critical Sobolev exponent. When Q is the unit ball it is easy to check 
that there exist radially symmetric positive solutions of 

Au + uP = 0, 

which have multiple blow up at the origin as the exponent p tends to pj\] (we do not assume 
Dirichlet boundary condition here). We discuss this result in section 4. In a recent paper 
[5], Del Pino, Dolbeault and Musso have proved that a similar result was also true for (1). 
These solutions which have multiple blow up at some points in O will be referred to as 
"bubble tree solutions". We are interested in the existence of these bubble tree solutions 
when SI is arbitrary. 

2 Statement of the result 

Let G denote Green's function for the Laplace operator with Dirichlet boundary condition 
on n and let H denote Robin's function, i.e. the regular part of Green's function. Namely 

for (y, z) G X ri. Observe that AyH = in SI x and G = on d{'^ x f2). 
Given m G N* and x := {xi, . . . , Xm) € S^"*, we define the mx m matrix 

M(x) := {mij)i<ij<m, 

whose entries are given by 

ma := H{xi,Xi) > and m^- := —G{xi,Xj) < 0, (2) 

if i 7^ j. Let p(x) be the least eigenvalue of M(x). We agree that /9(x) = — oo, if = Xj for 
some i ^ j. Finally, we define r(x) to be the unique eigenvector associated to p(x) whose 
coordinates are all positive and which is normalized so that its norm is equal to 1 (given 
the signs of the entries of M(x), it is easy to check that one can choose the eigenvector 
corresponding to the least eigenvalues to have coordinates greater than 0) . 

We define the open set 

P+ := {(A,x) G (M;)™ X ^T"}. 
Given /x e IR and (£i, . . . , G N'", we define 

by 

J^^iA, x) := A M(x) *A - cff ^ Af^ + C^^^ ^ log A^, 

i=l 1=1 
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where A = (Ai, . . . , A^) and where are two positive constants which only depend 

(i) 

on N and which will be defined in section 8. In the following, we denote some positive 
constant which only depends on A^. 

Finally the parameter e > is defined by 
Granted the above definitions, our result reads : 

Theorem 1 Assume that N >h and ji are fixed. Let (A, x) G he a nondegenerate 
critical point oj T^. Then, there exists £o > O'l^d for all e G (0,£:o) there exists Up a 

N-4 

solution of (1) with A := /x£Jv-2 and p = pn + £, such that 

m 
i=l 

(3) 

in the sense of measures, where the constant Cj^ is given by 

JV+2 

In other words, the sequence Up converges to (in any topology) away from the points 
Xi, as the parameter p tends to pN- Near each Xi the solution Up has multiple blow up in 
the sense that there exists c > (independent of p) , Xi^p G and parameters dij^p^i^ > 
such that ^ 

- < dij^p^^ < c, 

and 



lim 

P^PN 



Up{- + Xi,p)-{N{N -2)) 




= 0. 



Here 

and ro > is fixed small enough. Moreover, there is a relation between the parameters 
di,i,p,ij, and the parameters Aj since 

p^}^^^f'di,i,p,^ = {N{N-2))\. 

We briefly describe the plan of the paper. In section 3, we give some applications and 
some comments. In section 4, we recall some well known fact about radial solutions of 
Au + -u^ = when the exponent p is larger than the critical Sobolev exponent pN- In 
section 5 and 6 we give a new proof of existence of radial solutions. This proof is needed 
just because, for the proof of Theorem 1, we need some estimates which arc not available 
in [5]. Finally, the proof of the main result is the content of the sections 7 to 8. This 
proof is based on a gluing technic already used by Mazzeo and Pacard [10] in a different 
context. 
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3 Applications and comments 



Application 1 We consider the case where m = 1 and Q, = B \s the unit baU, we recover 
the result of Del Pino, Dolbeault and Musso [5]. Indeed, given ^ G N, we have 

•^^(^I'^i) = n _ i^i2w-2 -^'^N + eg) i log Ai 

It is clear that, provided the constant /x is chosen sufficiently large, this function admits 
two nondegenerate critical points which we denote by (Ai^i,0) and (Ai^2,0). Therefore, 
for any /x large enough, we find two distinct solutions of (1). 

Application 2 Now assume that is "close" to the unit ball. Then, a standard pertur- 
bation result shows that, for a given £ and provided ^ is sufficiently large, the function 
also admits two non degenerated critical points. This fact again guaranties the existence 
of two distinct solutions of (1). 



Application 3 We consider the case where m = 2, fi = 0, £i = £2- When CI = R^— S(0, 1) 
the functional J^o can be explicitely written as 

A^ A^ 

.^0(Al,A2,Xl,X2) = (1 _ |2)JV-2 + (1 _ I j2)jV-2 + '^N ^1 log(AiA2) 

-2Ai A2 ' ^ ^ 



It admits a critical point (A^, Ag, xj, Xg) where 

a; = a^ = 

and 



1/2 



2 2 2_ 



= (a*,0,...,0) = (-a*,0,...,0) 



where a* > 1 satisfies 



1 



+ 



(2a,)^-i (a2- 1)^-1 (a2 + 1)^-1' 

These explicit formula allow one to study (1) for /x = in a annular domain 0, = B{0, 1) — 
B{0,p), when p is close to 0. Indeed, we write z = {zi,z') G M x M''^"^ and using the 
symmetries, it is enough to look for solutions of (1) which only depend on zi and \z'\ and 
blow up at two points (which turn out to be close to dB{0,p) as p tends to pn). For 
this purpose, we study the functional J^q, reduced by the symmetries we impose. In a 
neighborhood of {p^~ Ai, p^~ A2, pa^,, pa^,) G M^, J^o can be expanded as 

2-N/ Af A2 . 



J^o(Ai,A2,s,t)= p'-''{ 



(1 - (l_(i/p)2)iV-2; 

_0A A ^ /^^-^ 

+Cg) £1 log(AiA2) + 0{Al + Ai) 
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And this functional admits a non degenerated critical point, provided p is sufficiently 
small. Applying the result of Theorem 1, we find solutions of (1) which have two bubble 
trees located near dB(0, p). When li = I2 = I such a result has been obtained by Felmer, 
Del Pino and Musso in [7] (see also [8] and [6]). 

Comments When m = 1, a necessary condition for the existence of critical points of J^fj^ 
is given by : /i is a sufficiently large positive number. Indeed, a nonexistence result for 
single peaked solutions of (1), when = 0, has been proved very recently by Rey et al in 
[2]. 

When m > 2 and /i = 0, if J-q admits a nondegenerate critical point, then Tfj_ also 
admits a nondegenerate critical point, provided // is small enough. This means that even 
for negative values of )U, we can construct solutions of (1). 

Also observe that in the case where = 0, if (A°,x'^) is a nondegenerate critical 
point of for (^1, . . . then (\/fcA°,x°) is a nondegenerate critical point of !Fq for 

{Ml, kim), where k eN. 

Finally, observe that our result parallels the corresponding result which has been ob- 
tained by Bahri, Li and Rey [1] for the subcritical case, i.e. when p < p^. In such case, 
only simple bubbles can be appeared, i.e. there are no bubble-towers (see [9]). 

4 Positive radial solutions of Aw + = in 

We recall some well known facts about positive radial solutions of 

/^u + vF = 0, (4) 

in M^. It is standard to look for radial positive solutions of (4) of the form 

2 

n(a;) = p-i ■«(— log |x|). (5) 

If we set t = —log I a; I, then v is a solution of an autonomous second order nonlinear 
ordinary differential equation : 

dlv - Up dtv -bpv + vf = 0, (6) 
where the constants ap and bp are given by 

4 J , ^ fl^T 



ap:=N-2 -, and bp := ^- . (7) 

p — 1 p — 1 \ p—i-J 

Observe that ap vanishes precisely when p = pN and bp vanishes when p = ;j^/z2- 
introduce the function 

Hp{x,y):=ly^-^-^x^ + ^^. (8) 
If V is a solution of (6), then 

dtHp {v, dtv) = ap {dtvf' . 
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In particular, this implies that dtHp {v,dtv) > when p > pn- 

There are two stationary solutions of (6), the first one is given by v = and the other 

1 

one is given hy v = bp~^ . We set 

1 

Cp '.= bp 

We claim that there exists a heteroclinic solutions of (6) when p > p^. This is the content 
of the following : 

Proposition 1 Assume p > pN- Then, there exists a unique solution Vp of (6) which is 
defined on R, satisfies 

lim Vp(t) = Cp, lim Vp(t) = 0, (9) 

t— >— oo t— >+oo 



and is normalized so that 



lim ep-i u„m = 1. (10) 



This solution satisfies H{vp, dfVp) < 0. 

Proof. We first prove that there exists a unique solution of (6) which is defined for t large 
enough and which satisfies (10). According to a classical result in the theory of nonlinear 

ordinary differential equations [4], it is enough to check that there exists a solution of the 
homogeneous problem associated to the linearized ordinary differential equation at = 0, 
which has the desired behavior as t tends to +oo. Now, the associated homogeneous 
problem reads 

d^v — Qp dfV — bpV = 0. (11) 
And clearly it has two independent solutions which are given by t — >■ e'''='=* where 

7+:=iV-2-^, and 7_ := ^. (12) 

p— I p — ^ 

Therefore, there exists a unique solution of (6) which is asymptotic to t — > e'"'-* as t tends 
to +00 and hence satisfies the second formula of (9). A priori this solution, which from 
now on is denoted by Vp, is only defined for t large enough, say t G (t, +oo). Observe that 
there also exists another solution of (6) which is asymptotic to t — > e'^+* as t tends to 

— DO. 

Since the function t — > Hp{vp, dtVp) is increasing and 

lim Hp {vp, dtVp) = 0. 

t — *4-cx) 

we conclude that Hp{vp,dtVp) < for any t G {t,+oo). Thus, Vp remains bounded inde- 
pendently of the value of t and hence can be extended to all M. Now, as t tend to — oo, 
there two possibilities : either Vp converges to a limit cycle or Vp converges to the constant 
Cp, the unique stationary point in region {{y, dtv) : Hp{y, dtv) < 0}. But dtHp{v, dfv) > 
if dtV 7^ 0. Hence, there are no limit cycle. We conclude that lim.t^_^Vp = Cp. This 
completes the proof of the result. ■ 
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In the next result, we show that the function dtVp vanishes at infinitely many points, 
provided p is close enough to pN- 

Proposition 2 Assume that p > pN and further assume that 

al - A{p - l)hp < 0. (13) 

Then the set of zeros of dtVp is given by two sequences (£i)i>i and (ii)i>i tending to — oo 
and satisfying 

U > U > ij+i > t^_^_i, 

Moreover, we have 

Vpiti) <Cp< Vp{ti). 

Proof. We linearize (6) at v = Cp. This yields the operator 

Lp = -apdt + {p- 1) bp, 
since c^~^ = bp. The characteristic roots of Lp are given by 

7± = ^ (flp ± ^ ^4(p - l)bp - aj). (14) 

These are imaginary valued since (13) is satisfied. It follows from standard theory for ordi- 
nary differential equations [4] that Vp — Cp is asymptotic to a solution of the homogeneous 
system associated to Lp. Hence there there exists c, d G IR and 7 > ^ such that 

Vp = Cp + c^{e"^+'+^) + 0{e^'), (15) 

as t tends to — oo, where 5R(-) is real part of a complex number. This immediately implies 
that dtVp has infinitely many zeros. The result of the proposition follows at once from this 
expansion. ■ 



We define 



We now derive an upper bound for the solution Vp which has been defined in Proposition 1. 
This upper bound follows from the more general result : 

Proposition 3 Assume that v is a solution of (6) such that Hp{y,dtv) < on {ti,t2). 
Then \v\ < dp on (^1,^2)- 

Proof. This follows at once from the fact that 

max{x > : 3 y G M Hp{x, y) < 0} = dp. 
together with the fact that we have assumed that Hp{v, dtv) < 0. ■ 
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Prom now on we assume that (13) is satisfied and we define the sequences 

£p,i = Vpiti), and rjp^i = Vp{ti), 

which correspond to the sequence of local minima and local maxima of the function Vp. 
Observe that we have the sequence {£p,i)i (resp. {rjp^i)i) is increasing (resp. decreasing) 
and converges to Cp 

< Sp^i < ep^2 < ■ ■ ■ < Cp < . . . < r/p,2 < Vp,i < dp- 

It will be convenient to agree that 

to = +00 and £pfi = 0. 

We now derive a precise expansion of the value of Ep^i as p tends to the critical exponent 
Pn. This result relies on the following more general result which gives the asymptotic of 
the first return map when p is close to pat- 

For p > Pn and rj G [£p^i,Cp], we consider the function Vp^r] which is a solution of (6) 
which is defined in (0, tp ,^) and satisfies 

Vp,r,iO) = rj and dtVp^rj{0) = dtVp^r){tp,r^) = 0. 

If tp ,^ = +00, we agree that the above equalities have to be understood as limits. We fur- 
ther assume that Vp^r) is strictly increasing on (0, tp_^) and strictly decreasing on iip^n,tp,ri)- 
Finally, we assume that 

Hp{vp^n, dtVp,r)) < 0, 

on {0,tp jj). In other words, is the first return time. Observe that, when p = pN the 
equation satisfied by Vp^r) is Hamiltonian hence we have 

pl™^(^'p,r?(0) - Vp^r,{tp,n)) = 0- 

We make this estimate more precise in the following : 

Proposition 4 There exists a hounded positive function Dn '■ [0, Cp^,,] — >■ M"*", which only 
depends on N such that 

lim ^P.^(O)' - ^ 

p^pn P — Pn 

uniformly with respect to rj. 

Proof. For the sake of simplicity in the notations, we drop the p, 77 indices. Since we have 
assumed that Hp{v,dtv) < in {0,t), we get 



Recall that 



\dtv\ <JbpV^ -^^'"^^^ < V^^- 



dtHp {v, dfv) = ap{dtvf'. (16) 
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Integrating this equality over {0,t) and using the fact that dfV > on {0,t), we get 

< Hp{v{t),0) - Hp{v{0),0) = ap f {dtvf dt 

Jo 



„i 

< ttp \ hp i dtv V dt 

(17) 

^ Jo 

Similarly, using an integration over (t, t) together with the fact that dtv < over this set, 
we also get 

< Hp{v{V), 0) - Hp{v{t), 0) < ^ Op ^pvit)". (18) 
Hence, we conclude that 

< H{v{t),Qt) - H{v{Q),Q) < ap ^Jbpv{ff. (19) 

Thanks to the previous Proposition we know that v{t) < dp and clearly 

lim ^ = (20) 

while bp and dp remain bounded as p tends to pN. This, together with (18) and (19), 
implies that 

^hm {Hp{v{0),(}) - Hp{v{t),0)) = ^hm {Hp{v{0),0) - Hp{v{t),0)) = (21) 

uniformly with respect to rj. As a consequence, we get using the expression of Hp the fact 
that 

lim (v(Of - v(tf) = lim (vq - v(t)) = 

uniformly with respect to rj, where vq > v{0) is the unique solution of Hp{vo,0) = 
Hp{v{0),0) which belongs to {cp,dp). 

This being understood, we write 



/* {dtvf dt = f J2Hp{v{s), dtv{s)) + bpv^- ^- vP+^ dtv ds 
Jo JO V J* + 1 

rit) I 2 

= / X '^Hp(x,dtv(x)) + bpx^ xP+^dx. 

Jv{0) V P+'^ 

Now, as p tends to pn, it follows from the previous discussion that the right hand side 
converges (uniformly with respect to rj) to 



n {N-2V N-2 2N 

EN{r]) := ^2Hp^{r],0) + ^ ^ ^ - -^x—^dx 

where fj > Cp^ satisfies Hpj^{r],0) = Hp^{fi,Q). Similarly, we have 

lim f~ {dtvf dt = ENiv) 



where the convergence is uniform with respect to rj. Moreover the function ij — > E]\f{rf) is 
bounded. Using these Hmits together with (16), which we integrate over (0, i), we conclude 
that there exists a constant E]\f{r]) := ^^~2^^ EN{r]) only depending on N such that 

p-^pn p — Pn 

uniformly with respect to rj. The result follows at once from these limits and the expression 
of Hp. m 

Looking at the previous proof, it should be clear that 
Proposition 5 As p tends to pjv, the functions 

'Vp,ri •= '^p,ri{' + tpjri) 

converge (uniformly on compacts) to Wpj^^r] the unique solution of 

d^w -bp^w + = 

with w{0) = fj and dtw{0) = where fj > Cpj^ satisfy Hpj^{rj,0) = Hpj^{fj,Q). Moreover, 
the convergence is uniform with respect to rj. 

Proof. This follows at once from Ascoli's theorem since Vp^.,^ and all its derivatives are 
uniformly bounded. ■ 

Observe that, in the previous Proposition, as r) tends to 0, the function Wpj^^r] converges 
(uniformly on compacts) to wq which is explicitly given by 

wo{t) := -) (cosht) 2 . 

Going back to the study of the function Vp, the result of Proposition 4 yields: 

Corollary 1 There exists a positive constant C^^ (in fact C^^ = Dn{0) given in Propo- 
sition 4), only depending on N, such that, for all i eN 

lim -A^=iC^^) (23) 

P^Pn p - pi^ 

Moreover, we have the explicit formula for 



(4) _ ( N{N-2) \ (AT -2)^ /■+~ dt 

\ A J 2{N-l)J-oo (cosht)^-2 ^ ^ 



In the next result, we estimate any solution of (6), near one of the points where it 

achieves a minimum, by comparing it to the solution of a linear problem. Indeed, we 
consider Wp to be the solution of the second order linear ordinary differential equation 

dfWp — Qp dfWp — bpWp = 

(25) 

wp{0) = 1, dtWpiO) = 0. 
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which is exphcitly given by 



where 7± have been defined in (12). The following Lemma shows that, close to 0, the 
solution Vp^ri of (6) with ■Up^^(O) = rj and 5ti^p,,j(0) = is well approximated hy rjWp. 

Lemma 1 For all k eN, there exists a positive constant > such that for all t eM. 

\d^{vp,r,-vyJp)\<CkV''wP (26) 

for p close enough to pN- 

Proof. Again wc drop the indices p, r/ to keep the notations simple. We view t; as a solution 
of a non homogeneous linear second order ordinary differential equation 

dfV — Op dtv — bpV = —v^ 

The variation of the constant formula yields 

v{t) = r) w{t) - e^+* /* e("^-2T+)^ e(-"^+^+)^ v{Cy dC ds, (27) 
Jo Jo 

This in particular implies that v{t) < r)w{t) for all t G R. 

When i > 0, we can therefore use the bounds 

v{t) <r]w{t) <cr]e'^+^ (28) 

in (27) to conclude that 

\v{t)-r)w{t)\ < 07^6^+* /■*eK-27+)^ e^-''r'+i+)i e^+P^ dQ ds < erf e^^+K (29) 

Jo Jo 

When t < 0, a similar analysis yields 

\v{t) -r]w{t)\ < cTfeP^-* for all t G (-oo, 0). (30) 

This completes the proof of the estimate of v. The estimates for the derivatives follow 
similarly. ■ 

The last result translates for the function 

2 

'^P,vi^) ■= \x\ %r?(- log 

and we obtain the estimate 



7+ 7- I |2-iv 



<Ckrf {\x\-P^- + \x\-P^+) \xp-, 



N-2' 

where the constant Cfc > only depends on k and N and remains bounded as p ^ pn- 

As a consequence, we have the following result which provides an expansion of t^ and 
ti as p tends to pN '■ 
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Corollary 2 As p tends to pjv, we have 

2i-l 



ti — — log £ 

' N -2 ^ 



+ 



loge 



< Ci, (31) 



N-2 

for some constant Cj > which only depends on N and i. We recall that e = p — pN- 

Proof. As t goes from to t^+i, the function Vp passes once through the value Cp. Hence 
there exists G such that = Cp. 

We first estimate — U+i- In view of the previous Proposition, this quantity can 

be estimated by 

- ti+1 = ^-^loge + 0{l) (32) 

Now, we claim that tj+i — remains uniformly bounded as p tends to pjy- Indeed, 
it follows from the remark after Proposition 5 that, as p converges to pN, the sequence 
of functions t — f (ti+i + t) converges on compacts to wo(t) = ^ ■^(■^~^) (cosh t) 2 . 
From this we conclude that it takes a finite time for wq to go from Cp^^ to dp^. Hence, 
provided p remains close to pN, the time it takes to Vp to go from Cp to Vp(ti-^-l) is bounded 
uniformly as p tends to pN- 

Therefore, we conclude that 

U+i-ti+^ = --^logs + 0{l). (33) 

Similarly, we obtain 

t,-ti+i = ^loge + 0{l). (34) 
In order to obtain the estimates as stated, just observe that 

and also that ti = 0{1). ■ 

Now, we compare solutions of (6) which have different boundary data. We keep the 
previous notations. We prove the following technical result : 

Lemma 2 For all cq > 1, there exists a positive constant c > only depending on N and 
Co such that 

^ (77 - r/) < \vp,r,{tp^r,) - Vp,fiitp^ii)\ <c{fi-ri) (35) 



and 



\tp,r, - tp,fj\ < c ^ (36) 



for all p close enough to pjv, provided (^ + ^fc^)£2 < r) < fj < {cq + \fc^) £2 
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Proof. We set v = Vp^jj and v = Vp^fj- To prove the result we write for the difference 
D := V — V 

d^D - GpdtD -bpD = ~fD 

where 

^ _vP-vP 

V — V 

It follows from the estimates of Lemma 1 that, for all p close enough to pN, 

\f\<c{vwpr-^ 

for some constant c which only depends on A'^ and cq. Now, as in the proof of Lemma 1, 
we use the variation of the constant formula to get 

D={fj-v)wp-Wp f e«^^ w-\s) r e-«-^ WpiQ /(C) 0(0 dC ds (37) 
Jo Jo 

We are interested in the range of validity of the two sided estimate 

^ifi-ri)wp<\D\<2{fj-ri)wp (38) 

Inserting this into (37), we get 

{fi-ri){l-c{riwpf-^)wp< \D\ < {fj - rj) {1 + c {rj Wpf'^) Wp 

Form which it follows that (38) is valid up to the time tp^rj where c{r]WpY~^ = 1/2. 
Therefore, we have 

\{fl-ri)<r] \D{ip,r,)\ + r? \dtD{ip,r,)\ <2{fj-ri) (39) 

at this point. Now, it should be clear that ip^ri — tp,ri is bounded independently of for p 
close to PN (since {v,dtv) remains bounded away from in this interval). Hence we also 
have 

^ (?7 - ^7) < ^ l(^^ - v)iip,r,)\ + V \dt{v - v)iip,r,)\ <cifj-ri) (40) 
for some constant c > 0. Standard result on dynamical systems imply that 

V\tp,fj-tp,rj\ <c{fl-1]). (41) 

Using (16) and (40), we have 

\Hp{v{tp^r,),0) - Hp{v{tp^r)),dtV{tp^r)))\ 

= \Hp{r,, 0) - Hp{fi, 0) - ap J^"'" [{dtvf - (dtvf) dt\ (42) 
< cri{fj — rj) 

since (^ + ^/CN{4))e^ < r] < fj < {cq + \/Cjv(4)) e^. Together with (41), we estimate 

\Hpivitp,r,), 0) - Hp{v{tp,fj), 0)1 < cr){fj - r?) (43) 
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which implies 

-Viv -V)< \vitp,ij) - v{tp^r,)\ < cri{fj - ri). (44) 
Prom Corollary 1 and results on system dynamic, there holds 

2^ £5 < v{tp^^) < V%^f,) < 2co £5 . (45) 

Using similar arguments on {ip,rptp,r))i we get 

^ \{t,p,f, - tp,fi) - {tp,n - ip,r^)\ <c{fj-r]) 

and also that ^ 

-{fj -r))< \v{tp^i^) - v{tp^^)\ < c{fi - rj). 

The result follows at once from these estimates. ■ 



5 Linear results 



We keep the notations in the previous section. For the sake of simplicity in the notations, 
we drop the indices p and rj. We consider w to be the solution of 



d^w — Qp dtw — bpW + pv 



w = e 



-■2t , 



(46) 



in (0,t) with boundary conditions w{t) = dtw{t) = 0. We are interested in the behavior 
of w as f] tends to 0. This is the contain of the following result. 

Lemma 3 Assume that N > 5. Let cq > 1 and dg > be fixed. Assume that rj G 
(^e2jCo£2). Then, there exist c > and eq > such that for all p G {pn,Pn + £o) we 
have 



wit)- 



^Pp,ri 



r}{N - 2)2 

in {—do, do), where the constant (3p^r) is given by 



f^P,r] := i'p,r,(«) e~^^ ds. 



Moreover, we have 



P^PN ^" \ 4 



-2s 



(cosh s 



ds := C 



(5) 
N ■ 



(47) 



(48) 



(49) 



Proof. As usual we drop the p, rj indices. We use the fact that wi := rj ^ dfV is an explicit 
solution of the homogeneous problem 



dfWi — Up dtwi — bpWi + pv^ ^ wi = 0. 

This yields a representation formula for w, at least when t G {t,t). 

ft ri 
w{t)=wi{t) e^p'wi^is) / e-''^'^ wi{C) e~^^ v{C) dC ds. 

Jt Js 



(50) 
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Observe that the result of Lemma 1 yields 

ie^-(*-*) <^;(t) <ce^-(*-*) 

for all t G {t, t) and also 

c 

for all i G {t + l,t — 1). Using this, we get the estimate 

\w{t)\ + \dtw{t)\ < ce-2*eT-(*-*) (51) 

for some constant c > 0. This estimate is valid for all t G {t + l,t), however, enlarging 
the value of c if this is necessary, we can assume that this estimate holds for t G (? — l,t). 
The solution w extends to (0,t). 

Again, we use the fact that wi = r]~^ dtv and 

t-i 



W2 (t) := wi (t) e""' (g) ds 



which is defined for t G (l,t — 1), are solutions of the homogeneous problem (50). Hence 
we can decompose 

w = aiwi + a2W2 + w 

where w is defined by 

w{t):=wi{t) fe'^^'w^^is) e-""^ wi{C) e-^^ v{C,) dC ds. 
Jo Jo 

As above, the result of Lemma 1 yields 

-rje''+^ <v(t) <crie^+^ (52) 
c 

for alH G (0, i) and also 

^r/e^+* < \dtv{t)\ < cr?e^+* (53) 

for alH G (1, i — 1). Using this, we get the estimate 

\w{t)\ + \dtw{t)\ < cr?e-2*eT+* (54) 

for some constant c > 0. This estimate is valid for all i G (0,t — 1). 
Since 

- <r]\wi{t-l)\ 
c 

it follows at once from (51) and (54) that we can estimate the parameter ai by 

|«i| < C77e~^*. 
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In order to estimate the parameter we multiply the equation (46) by w\ and integrate 
by parts. Using the fact that wx is a solution of (50) we obtain 



\w\dtw — w dtW\ — apWW]\-^ = —2ap / dtWiwdt + r/ / e vdtvdt. 

Jo Jo 

Since wi = at t = and w = dtw = at t = t, this simplifies into 

w{0)dtwi{0)= -2ap [~dtwiwdt + rj^^ [~ e""^^ v dtv dt 
Jo Jo 

Prom (52) and (53), it follows 

\w2{t)\ + \dtW2{t)\ < ce-T+* 

for all t G (1, t — 1). Enlarging the value of c if this is necessary, we can assume that this 
estimate holds for t € (0,t). 

Collecting these estimates, we get 



/ 

Jo 



1 2 f— 1 2 

wdtwidt = a20{loge) +0{e~^+ir:2')^ and J_ wdtWi = 0{e'~^'^N^). 



To calculate ^2(0), we see that 

e''p\pvP-\t) - bp) 



1 d f e"f* 



dtv{t)dt \dfv{t) 



Hence we get 



W2{t) 



wi{t) y* \'''''w^^{s)ds + wi{t) J\''p'w^^{s)di 

It-i * 
wi{t) / e'''''wl^{s)ds + wi{t) 

-•i e^p'ibp-pvP-^is)) 



wi{s)dtwi{s) 



-Wl 



(dtm(s)y 



-ds 



for all t G (—do, do)- In particular 



W2{0) = 



since wi(0) = 0. Consequently, we obtain the estimate 



a2 



Jo 



ds + 0{ei). 



{wi{s)f 



It remains to estimate W2 in the neighborhood of 0. We first estimate j t^ttt- We 

decompose 
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It follows from (53) that 



J 

J12 



and 



1 < e^f^ < 1 + ce log-, 

£ 



for all s G (0,f). Using the result of Lemma 1, we obtain 

N -2 N -2 2 

dtv{t) = — - — ?7sinh( — - — + 0(eJv+2 )) 



(55) 



for all t E Ii. Therefore, we deduce 



On the other hand, using again Lemma 1, we have 

wiit) 



COSh(^) \ 2 



N -2 . ^ N -2 ^ P-i. 
_^smh(^— 1)+0(£— , 



dtwiit) = ^ ^ ' cosh(^— t) + 0(£— , 

for alH G (— do, f^o)- Now direct calculations lead to 

4 



(56) 
(57) 



W2{t) 



(N-2)t 2 

e" 2 +0{eN+2] 



{N - 2)2 

for alH G (-do, -c^o)- This proves (47). 

Finally, in order to obtain (49), it is enough to observe that v{t+-) converges (uniformly 
on compacts) to wq. This completes the proof of the result. ■ 



Using similar arguments (and the notations of the previous Proposition) , one can show 



that 



if AT > 6, and 



w{t) 



4/3 



ri{N - 2)2 



-{N-2)t/2 



1 I 3 

<cj3e 2 + iv+2 



< Cl3e 2 + iV+2 iv2_4 



r]{N - 2)2 

if iV = 5, for all t G (^^ log \ - do, log I + do). 
In the following, /3p^^ will be expanded. 

Lemma 4 Under the above assumptions, let co > 1 be given. Assume 



(58) 
(59) 



Then, 



where 



1 1 1 

— £2 < T] < CqS^. 
CO 



J^^og^ + cP+0{e^ log 



(60) 



"iV-2 



^ log2 + ^logA^(iV-2). 
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Proof. We recall 

t= / = (61) 

Jv j2Hjv,dtv) + br,v^-^ 

We divide 



:= hUhUh 



We estimate 



/ , = ={1+0{£n)) 

Jh ^2Hp{v,dtv)+hpV^-^ Jh ^bp{v^ - 772) 

= (1 + 0(6^)) ^ _ ^ Argch{e^ /ri) 

f dv 2_ f dv 
[ ^ = / ^ + 0(e'^) 

Jh ^2H,{v,dtv) + b,V^-'-^ Jh ^ ^N^Yy2 _ N^^^, 

2-N 



Recall w^it) = (M^^) V (cosht)^. We deduce 



where WQ{t) = e 2jv . Hence, the desired results yield. ■ 
We set v{-) = v{t + ■). We have the following result 

Lemma 5 Given cq > 1, assume ■u(O) G {dp — CQ£,dp). Then, there exists the constant c 
independent of p such that 

(JV-3/2)|t| 

\v{t) - wo{t)\ + \dtv{t) - dtWo{t)\ < cee 2 (62) 

for all t G {—i,t — t), provided p close to pN- 

Proof. We write for the difference D := v — wq so that 

d^D-bp^D = -fD + g (63) 

where 

f ■■= — and g = Qpdtv + w^'^ - w^. 

V — Wo 

Clearly, there exist some positive constants K and c independent of p such that 
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for allt e {K,t-t)U {-t, -K) and 

\9{t)\ <ce 

for all t e i-t,t- t). Recall 

\v{0) - wo{0)\ + \dtv{0) - dtwo{0)\ < COS. 

Hence, the desired result follows from the standard ordinary differential equation theory. 

■ 

As a consequence, we obtain immediately 

(7) 

Corollary 3 There exists a positive constant C)^ (only depending on N), such that 

where := '^^ . In particular, 

ep,e={isCi^^y/' + OieW). (64) 
We keep the notations introduced in section 4 and we define for all ^ G N 

f3p,t:= f^~\l{t)e-^^'-'-^Ut. (65) 

Thanks to Lemma 3 to 5 and Corollary 3, we conclude 

Corollary 4 There exists a positive constant C]y (only depending on N) such that 

2 / /si 2 1 iV-17/4 \ 

= {is)—^ (cP + 0(8 N log - + e^^)j (66) 
where := {ci^^)^^CPe-^^'^\ 

6 Radial solutions of Au + Au + = in the unit ball 

In this section we recover part of the result of Del Pino, Dolbeault and Musso concerning 
the existence of solutions of (1) in the unit ball. In doing so our aim is to derive precise 
estimates for these solutions which will be needed in the forthcoming construction. 

We begin with the definition of weighted spaces in cylindrical coordinates. These 
spaces are at the heart of our construction. 

Definition 1 Given (5 G M and — oo < ti < t2 < +og, the space C^{{ti,t2) x S^~^) 
is defined to be the set of continuous functions w G Cig^{{ti,t2) x S^~^) for which the 
following norm is finite : 

lkllcO((ti,t2)x5Jv-i) := l|e"'^^ w||L°°((ti,t2)x5^-i)- (67) 
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We would like to prove the existence of radial solutions of 

Au + Xu+\u\P~^u = in 5(0,1) (68) 

Using (5), we reduce the study of (68) to study of the nonlinear second order ordinary 
differential equation 

d^v - ap dtv -bpv + \v\P''^v + A e"^* t; = 0, (69) 

in (0, +oo). We keep the notations introduced in section 4 and we consider the linear 
operator 

Lp,ri ■.= dt - apdt-bp+p v^'^ . (70) 

We state, without a proof a result which will be proven in a more general context in the 
next section. 

Proposition 6 Assume that 5 G (— —^^^) is fixed. Then, there existed > 0, r/o > 
and c > such that, for all e G (0, eo); for all rj G (0, rjo) and for all f G Cg{{0,tpjj)), there 
exists a unique solution w G C^((0,tp,j)) of 

Lp,r)W = f, (71) 

in (0, tp^^) which satisfies 

w{tp^^) = dtw{tp^^) = 0, (72) 

with p = pn + £■ Furthermore, 

ll^llco <c||/|lco- (73) 

When tp^ri < +00, the existence and uniqueness of the solution of (71) is straightforward 
but the uniform estimate (73) requires some work. When tp ,^ = +oo, the boundary data 
(72) have to be understood as limits as tp ,^ = +oo. 

The next result will allow us to recover (part of) the result of Del Pino, Dolbeault and 
Musso [5] : 

Proposition 7 Assume that i £ N is fixed and that N > 5. Then, there exists eo > 
such that for all G M, for all ^ gM. and for all e G (0, eo), problem ( 68) with p = pN + £ 

N-4 

and X = lie ^-2 admits a solution which can be written in the form 



N-2 
2 



up,,,^ix) = iN{N-2))^ [yT^WJ ^'^^^ ^^^^ 
where o(l) converges uniformly to on B{0, 1) as £ tends to and where 

ej := dj {e'"^)^ 

for some parameters dj which are bounded from below and from above by some positive 
constant independent of e. Moreover we have the following expansion 

1 



= (fe)2 



6 2+ -I e 2 x ^-^ ; e 2 

2 (N-2r^ 



+0{eiri) 

2 
71 



(75) 



in B(0, 2r£) — 5(0, re/2), where := e^'^-i . Furthermore, Up^\^^ is positive in S(0, 2r£). 
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Proof. The proof is decomposed in several steps. We give the prove in the case where 
A'^ > 6 since, when N = 5, the proof is similar with straightforward changes. Given ^ G M 
(which will be fixed later on) we define 

T2i = U-i -te + ^y T2i-i = -t^ + ^ 

forO < i < £ and 

T2e = +oo, T2e-i = ti-t( + (, To = 0, 

For all < i < ^ — 2, we define Vp^i to be the solution of (6) in [T2i, T2i+2] with boundary 
conditions 

Vp,iiT2i+2) = £p,i-l-i + 0!i, dtVp^i(T2i+2) = 

and we define Vp^£-i to be the solution of (6) in [T2i-2, +oo) with boundary conditions 
Vp,e-i(T2e-2) = dtVp^e-i{T2i-2) = 0, 

for some parameters G M (which are assumed to be small). Here the parameters £p^i 
are the one which have been introduced in section 4. 

For any < i < £ — 1, we define the function 

Wi{t):=Vp^i{t + ti)+Wi{t) (76) 

for on the interval [T2i,T2i-\-2], for some parameters G M and some functions Wi G 
C^{[T2i,T2i+2])- We agree that t^-i = and ae-i = 0. 

Granted the above definitions, our strategy is the following : In Step 1 and 2, we 
look Wi solutions of (69) on each interval [T2i,T2i-\-2]- Moreover, Wi are positive if i > 1. 
In Step 3, we choose the parameters (ckq, . . . , 0^-2) and (to, • • • , ^^-2) so that the Cauchy 
data of Wi and of Wi-i coincide at T2i. Gathering the functions Wi together, we obtain 
a solution of (69) which still depends on ^. 

Step 1. For each 1 < i < £ — 1, we look for a solution of (69) in [T2i, T2i_|_2]. Recall that 
£ = p — Pn- We now assume that 

as £ tends to 0. We define the operator 

Lp^i = 5| -Updt- bp + pvP-\- + ti). (78) 
With these notations, the equation we need to solve reads 

Lp,iWi = -A e"^*(fp,j(- + ti) + Wi) - Qi{wi) (79) 

where we have defined 

Qi{wi) := \vp^i{- + ti) + Wi\P'''^{vp^i{- + ti) + Wi) - v^i{- + ti) -pv^'i^i- + ti) Wi. 
We fix the weight parameter 6 G (— "^^f^, consider the set of functions 

S^,i = {we Cl{{T2i,T2i+2)) : \\w\\co < KAe-('^+2)^2*+i} , 
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where the constant k> will be fixed later on. 

d+6 N-4 

Given w G £K,i, it follows from (31) that \w\ < ckXe^-^ . Recall that A = /x£Jv-2 and 
S > hence, we obtain 

I I 2+6+N 
\w\<CK£N-2 <^ Up -I- i-) 

on {T2i, T2j+2). Therefore, we are allowed to use Taylor's expansion \{l + ty — 1 —pt\ < ct^ 
for t close enough to 0, to estimate 



mw)\<cv;/{-+t. 



w 



Using this, we obtain 



Next, we estimate 



IIQiHIIco <cKf A£^^^^^^^e-(^+2)W. (80) 



since 7+ — 5 — 2 > and 7_ — 5 — 2 < provided £ is close enough to 0. With similar 
arguments, we get 

WXe^"^^ w\\co < Ae"2T2,||^||^^ < ^2 ^^^^ g-(5+2)T2i+i ^ ^^32) 

s s 

Combining (80) to (82), we have obtained 

\\Qi{w) + Xe-^^ ivp4- + ti)+w)\\co < cAe-('^+2)^2i+i (^f e ''^'Iv-"''"'' + 1 + k A e 1^ ) , (83) 

which holds for all w € iS^.i- 

Given w G we apply the result of Proposition 6 which provides a solution of 

Lp,i w = -Qi{w) - A e"^* {vp{- + U) + w) 

with w{T2i+2) = dtw{T2i+2) = 0. Thanks to (83), we also have the estimate 

MM IX 91T . „ (2+^+JV)(p-l) 4 ^ 

llu-llc" < A e(-*"^'^2i+i g (/.P £ +1 + K\e^). 

S 

for some constant c > which does not depend on w, nor on k nor on e provided this later 
is chosen small enough. This estimate being understood, we choose the constant «; > so 
that 

(2+a+jv)(p-i) 4 . 

c(K;^e Jv-2 + 1 + «;A£Jv-2) < AC, 

for all £ close enough to 0, say £ G (0, £o)- 

To summarize, using the above analysis, we can define the mapping 

by %{w) = w. Thanks to the above choice of k, the mapping % is well defined. Observe 
that this mapping is clearly continuous and compact so that one can refer to Schauder's 
fixed point Theorem to obtain the fixed point of %. We have proved the : 
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Lemma 6 Assume that ai and ti satisfy (77). Then, there exists Wi a positive solu- 
tion of (69) in (T2i,T2i+2) with boundary conditions Wi{T2i+2) = Vp,i{T2i+2 + ti) and 
dtWi{T2i-^-2) = dtVp^i{T2iJf-2 + ti)- In addition, we have the estimates 

\\Wi - vp,i{- + ti)\\co < c Ae(-'5-2)^2^+i 

o 

where the constant c is independent of e and of the parameters ai, ti and ^. 

Observe that the solution we have obtained is unique and depends continuously on the 
parameters a^, ti and ^ since it is the unique solution of an ordinary differential equation. 
This fact is even true when i = l—l even though the solution is defined on a half line. 

Step 2. We now look for a solution of (69) which is defined on {Tq,T2). We decompose 

Wo{t) = Vpfl{t + to) + w{t) + w{t), 

where w is the solution of 



Lpfi w = -Xe ^*Wp,o(t + to) 



in (To, T2) with boundary data w(T2) = dtw{T2) = 0. The operator Lp^Q is the one which 
has been defined in (78). With this in mind, it remains to find a w solution of 

LpfiW = — Ae~^*(uJ + w) — Qo{w) 

in (To, T2) with boundary data w{T2) = dtw{T2) = 0, where 

Qoim.) ■= \vp,o{- + to) +w + w\^^\vpfi{- + to) +w + w) -v^^q{- +to) -pVp^^{- + to){w + w). 
It will be convenient to define 



mm 



2iV + ^ - 6 {N + S-2)p \ 
N-2 ' J ■ 



Observe that we have q > ^ since we have assumed that N > 5 and S G (— — "^^2^)- 
This time we consider the following set of functions 

£^,0 = {we Cl{{To,T2)) : \\w\\co < Ke'^} , 
where the constant k> will be fixed later on. It is clear that 

||Ae-%,o||co < cAe(-'^-2)^i < cXe^ 
Using the result of Proposition 6, we get 

_2t i. N-2+S N-2+S 



\w\\(;0<c\\Xe Vpfl\\(;0<Cll£ ^^"2 <C£ ^^-2 
5 5 



so that 

As in Step 1, we have 



„. 2N-6+S 

\Xe uJ||co <c£ N-2 . (84) 

5 



|Ae~^*u;||co < X\\w\\co < ck£^+^. (85) 

5 5 
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for all w G 8^,0- For £ small enough, we have 1 < p < 2. Thus, for all S2 € K and all 
si > 0, we can write 



<C\S2\P. 



for some constant c > 0. Consequently, we can estimate for all w G £"^,0 



IIQoHIIcO < C SUPt^^T^^T^) e \w{t) + w{t)\P 



p(N-2+S) 



(86) 



Using the result of Proposition 6, we get a solution w of 

LpflW = — Ae""^* {w + w) — Qo{w) 
with w{T2) = dtw{T2) = 0. Collecting (84), (85) and (86) we get the estimate 

<e'ic{l + KFe^-^^'i + Ke^). (87) 
We choose the constant n so that 

for all £ is close to 0, say e G (0, eo)- 

As Step 1, we can define the mapping 

by Tq{w) := w. Clearly, Tq is well defined and is continuous and compact, so that one can 
again refer to Schauder's fixed point Theorem to obtain the fixed point of Tq. We have 
proved the : 

Lemma 7 Given ao and to satisfying (77), there exists a solution Wq of (69) in (To,T2) 
with boundary conditions Wq(T2) = Vpfi{T2 + to) and dtWo{T2) = dtVp^o{T2 + ^o)- In 
addition, we have the estimates 

\\Wo - Vpfl{- + to) - w\\co < ce^ 

o 

for some q > ^ and 

N-2+6 

\\w\\ro < ce ^-2 

II IIL,^ — 

Again this solution is unique and depends continuously on the parameters ao, to and ^. 

Step 3. We now explain how to choose the parameters (ao, • • • ,0^-2) and (to, ■ ■ ■ ,i£-2) 
so that the Cauchy data of Wj and Wi_i coincide at T2i, for 1 < i < £ — 1. To this aim, 
we argue inductively, starting by matching the Cauchy data of W^_2 and VF^-i- 

This amounts to find a£_2 and t^_2 so that 

W^_2(r2^_2) = W^-i(r2^-2), and dtWe-2{T2e-2) = dtWe-i{T2e-2), 
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In other words, we need to find Q!^_2 and t^_2 so that 

Vp/-2{T2e-2 + te-2) = Wi-i{T2t-2) dtVp^i-2{T2t-2+ti-2) = ^t^,?-! (72^-2) 

It follows from Lemma 1 and Lemma 6 that 

S+N+4e-(i 

W^-l(T2£-2) = £p,l+We-l{T2i-2) = £p,l+0{£ iv-2 ) 

S+N+4e-6 

dtWe-iiT2e-2) = dtwe-i{T2e-2) = 0{e ^-^ ). 
and we also have 

Wl_2{T2t-2) = (ep,l+«^-2)cOsh(^t^_2)+i^^-2 

dtWi_2{T2i-2) = {ep,i + a^-2)^ sinh(^i^_2) + G,_2, 



(89) 



where the continuous functions (Ff^2, Ge-2) depend on a^_2 and t£-2 and satisfy |-F^-2| + 

P 

\Gi-2\ = 0(e2). The system (88) is therefore equivalent to 

ti^2 = £^^Fl-2 and ai_2 = G;_2 (90) 

where the continuous functions (-F)_2,G/_2) depend on a£^2 and t^_2 and satisfy |F£-2| + 
|G£_2| = 0{e2) (here we have used the fact that ^^^^2~^ > 2 Provided e is close to 0). 

Recall that p = pn + £■ Given 7 G (1, |) we define 

B := {{ae-2,ti-2) e : a'-2 + £i'-2 < e''} ■ 

In view of (90), the mapping 

H :B — > B 

defined by H{ae-2, te-2) '■= (e"^ F^_2, Ge-2) is well defined and it follows from Browder's 
fixed point theorem that (90) admits a solution. In addition, applying Lemma 1 and 
Lemma 2, we get 

Vp,i-2{T2e-4 + te-2) = £p,2 + 0{e^), dtVp,e-2iT2e-A + te-2) = 0{e^. (91) 

Arguing inductively, we construct a function Wp,A,5) solution of (69), which depends on 
A and on ^, and which satisfies 

VpxdT2) = ep,i-i + O(e^), dtVyXi^T'2) = C(ei). (92) 

In view of (58) and thanks to Lemma 1, Lemma 7, Corollary 3 and Corollary 4, the 
following expansion holds 



Cj^^ (N-2)a-t) yCj^^ (2-iV)(g-f) AfJ,CP£^^-2 {N-6)i (2-N)t 
6 2 -|- — g 2 — , 6 2 g 2 



(iV-2)VcL'^ 



(93) 



(r. 1 I 3 1 (Af~2)t 1 2 -I JV-17/4 (2-N)t 

^ggdt ^ g-2 + iV+2 + (£2e 2 )P + £2 (£lV log i + £ N-2 )e 2 

for all t G (log log ^). This completes the proof of the result. 
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7 The linear analysis. 

Assume that is a regular bounded open subset of and S := {ai, . . . , am} is a finite 
set of points of CI. We choose ro > in such a way that the closed balls B{ai,2ro), for 
i = 1, . . . ,m are disjoint and included in Q. For all r G (0, ro), we define 

m 

^int,r '■— and ^ext,r '■= ^ — ^int,r- 

1=1 

We define the weighted spaces : 

Definition 2 Given S M, the space C^{Q — S) is defined to be the set of continuous 
functions w G Ci^JVi — S) for which the following norm is finite : 

m 

ll^llcO(Q-S) — + H Wr'^'wiaj + ■)\\L^{B{Q,2ro)-{0])- (94) 

j=l 

Given r G (0, ro) we define the space C^{^ext,r) to be the space of restrictions of functions 
of C^(fi — S) to ^ext.r- This space is endowed with the induced norm. 

In this section, we study the linearization of the nonlinear operator (68) about the 
radial function 

2 

Ue{x):=\x\ p-^ Vs{—\og\x\) 

where t;^ := Vp \^^ and Vp^\^^ is the solution of (69) defined in Step 3 of the proof of 
Proposition 7. This operator is defined by 

L, := A + A + ptxf-^ 

2 

Recall re = e'^'-'-^ . We can write any function v defined in the punctured ball .6(0, r^) — {0} 
as 

2 

'u(x) = |x| p-^ w{—log\x\,6), 
so that the study of reduces to the study of the linear operator 

Ce ■.= df- ap dt-bp + AgN-i + p vP'^ + A e'^* (95) 

on the half cylinder [Be, +oo) x where A^at-i denotes the Laplace-Beltrami operator 

on the sphere 5"^^^ and Be = — logr^. 

We denote by (e^, \j) the set of eigendata of A^jv-i 

Agjv-iCj = -\j ej. 

We also assume that the eigenvalues are counted with multiphcity, that \j < Aj+i and 
that the ej are normalized by 

/ e? da; = 1. 

JSN-I J 

We now prove some uniform estimates for a right inverse for the operator Ce- 
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Proposition 8 Assume that 6 G {—^y^,—^) is fixed. Then, there exists po G {pN, +00) 
such that, if p € {pn,Po), then, for all f G Cg{[Bs, +00) X S^~^), there exists a unique 
solution w G C'^{[B^, +00) x S^'^) of 

C,w = f (96) 

in [i?£,+oo) X S^~^ which satisfies 

w{Be, 9) G Span{ej : j = 0,..., N}. 

Furthermore, 

lkllc?<c||/||co (97) 



for some constant which does not depend on e. 



Proof. The proof is divided in three parts. In the first part we explain how to solve the 
equation (96) when the function / docs not have any component on ej for j = 0, . . . , A'^ in 
its eigenfunction decomposition. Next, in the second part, wc obtain a uniform estimate 
for the solution already obtained. Finally, in the last part, we explain how to solve (96) 
when the eigenfunction decomposition of / has components on eo, . . . , e^v- 

Step 1 For the time being, we assume that the eigenfunction decomposition of the function 
/ is given by 

f{t,e)= Yl fit) (98) 

j>N+l 



Observe that, as p tends to pN we have 



lim a„ = 0, lim br, 

P^PN P^PN 



N-2 



and 



hm^suppv;^,^^ = 



these limits being independent of the solution Vp^\^^. 

Now the eigenfunction decomposition of the Laplace-Betrami operator on 5^"^ induces 
a decomposition of the operator into the sequence of operators 



Lj ■.= df - Qp dt-bp - Xj + p v^~^ + A e~^*. 

Using these above limits together with the fact that Xj > 2iV for j > N + 1, we 
conclude that, for j > AT + 1 the potential is negative provided p is close enough to pN- 
In particular, this implies that it is possible to solve 

CeW = f 

on any (S^jS*) x S'^~^, with w = as boundary data (observe that the operator C/. is 
not self adjoint but is conjugate to a self adjoint operator and we have just seen that this 
former is injective, when restricted to the set of functions spanned by ej, for j > N + 1). 
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It remains to prove that there exists a constant c > which does not depend on S, 
nor on p such that 

sup|e-'^*ti;| < c sup|e-'^*/|. (99) 

Then, the existence of the solution on all {B^, +00) x S^~^ as well as the relevant estimate 
will follow by passing to the limit S — > +00. To keep the proof short and since anyway 
our aim is to pass to the limit as S tends to 00, it is enough to prove that (96) holds for 
all S chosen large enough so that sup(5 .,.go^ Vp < e. 

Step 2 The proof of (99) is by contradiction. If it were false for all choice of po and 
S, there would exist a sequence {pn)n tending to pN, a sequence of functions (/„) and a 
sequence of reals (S'n)n and a sequence {wn)n of solutions of (96) such that 

||/„|Uo = 1 and lim An := supe~^*|u'„| = +00. (100) 

■5 n— >+cx) 

We denote i?„ = i?e„ where := pn — Pn- Obviously, there exists a point {tn,0n) £ 
{Bn, Sn) X S'^^^ where the above supremum is achieved, namely An = e"''*" \wn{tn,0n)\- 
Observe that elliptic estimates imply that 

supe-''*|Vwn| < c{l + An) (101) 

and this in turn implies that the sequences (t„ — -B„)„ and {Sn — tn)n remain bounded 
away from 0. 

We define t„ > Bn to be the nearest local maximal point of the function Vp^ (t) to the 
point tn- We distinguish several cases according to the behavior of the sequence {tn)n- 

Case 1. Assume that the sequence {tn — in)n is bounded. In this case, we define the 
function Wn by 

Wn{t, 0) = ^ e--'*" Wn{t + tn, 9). 

Observe that the sequence of functions {vp„{- + tn))n converges on compact to t ^ {N{N — 
2))^^ (cosh i)~2~ (see Proposition 5). Up to a subsequence, we can assume that the 
sequence {tn — in)n converges to too- Moreover, we can assume that the sequence {'Wn)n 
converges on compacts to a nontrivial solution of 

-.9 . {N - 2f N{N + 2) , , , , , , 

dfwoo + Agiv-i Woo - — w'oo + — (cosh t)-^ Woo = 0. (102) 

Moreover, Woo is bounded by a constant times e^*. The fact that Woo is not identically equal 
to follows from the fact that \wn{tn — in,dn)\ = e'^^*""*") and hence remains bounded 
away from 0. 

We consider the eigenfunction decomposition of Woo 

oo 

Woc,= ^o^r 

j=N+l 

At — OO the function aj is either blowing up like t — > e~'^^* or decaying like t — > e'^^*, 
where 

r {N - 2)2 
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The choice of S E {—^^, — ^) imphes that —5 < -yj for all j > A'" + 1. Hence aj decays 
exponentially at — oo. Multiplying the equation (102) by ej and integrating by parts 
over M (all integrations are justified because aj decays exponentially at both =boo), we get 

oo ^ ^ J —oo 

J —oo 



< 



Since j > A?^ + 1, we have Aj > 2iV, and hence we conclude that aj = 0. Hence, i^oo = 0, 
a contradiction. 

Case 2. Assume that the sequence (tn — ^n)? the sequence {tfi — I^ji)n and the sequence 
{Sn — tn)n ai"© all Unbounded. In this case, we define the function Wn by 

Wnit, 9) = ^ e--^*" Wn{t + i„, 9). 

Observe that this time the sequence of functions (wp„(- + tn))n converge to on compacts. 
Up to a subsequence, we can assume that the sequence {wn)n converges on compacts to 
i^oo a nontrivial solution of 

(A^ - 2)^ 

dtWoo + A^jv-i Woo ^ Woo = 0. 

Moreover, Woo is bounded by a constant times e^^ . 

Again, we consider the eigenfunction decomposition of f^oo 

oo 

Woo= aj Cj 

j=N+l 

and we see that aj is a linear combination of t — > e~'^i* and t — > e^^^. The choice of 
d G (— , — -f-) implies that 5 > —jj for all j > A'' + 1. Hence aj cannot be bounded by 
e''* unless it is identically 0. We conclude that aj = 0. Hence, Woo = 0, a contradiction. 

Case 3. Assume that the sequence {tn — Bn)n is bounded (resp. that the sequence 
{Sn — tn)n is bounded) and that the sequence (i„ — in) is unbounded. This case can be 
treated as in case 2. The only difference is that this time Woo is defined on \t^, +oo) x S^~^ 
(resp. on (-oo,foo] x S^~^) and is equal to on {t^} x S^~^ (resp. on {too} x S^~^). 
We omit the details. 

Since we have reached a contradiction in each case, the proof of the claim is complete. 
We can now pass to the limit as S tends to +oo and complete the proof of the result in the 
case where the eigenfunction decomposition of / does not involve any Cj for j = 0, . . . , N. 

Step 3. Now we consider the case where the function / is coUinear to Cj, namely 

f{t,e)^f-'{f:)e,{e) 

for some < j < N. We extend the function / to be equal to when t < and we 
define the function Vp which is equal to Vp for t > and is equal to fov t < B^. We 
consider the equation 

dfa^ - ap dta^ - (Aj + bp) a^ +pvP-^ + Ae^^* x = f (103) 
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in R. Here x is a cutoff function identically equal to 1 on (^£,+00) and equal to on 
{—oo,B^ — 1). Observe that 

\m\<\\f\\coe''. 

For p close enough to pN, S is not an indicial root of the operator Ce and it follows from 
Cauchy's theorem that there exists a unique solution of (103) which is bounded by a 
constant times e'^* at +00. A priori this solution is only defined for t large enough but 
is can be extended to all R easily. Furthermore, it follows from the construction of this 
solution that 

sup e-^*y\ <csupe-^*\f\ 

(T,+oo) R 

provided T is large enough. This solution satisfies 

d^a^ - ap dta^ - {Xj + bp) = (104) 

for i < — 1 and, since S G {—^^^, ~t)' ^ven if blows up at —00, it blows up at a 
slower rate than t ^ e^^, provided p is chosen close enough to pN. 

We claim that there exists a constant c > such that 

supe"''* I a'' I < c supe"''* |/^ | 

provided p is close enough to pN- As before, we argue by contradiction. Assume that the 
claim is not true. Then there would exist a sequence {pn)n tending to pjv, a sequence of 
functions {f^)n and a sequence of solutions (a^)„ of (104) such that 

supe"'^* \ fl\ = 1 and An := supe~^* |a^| 

tends to +00. The previous remarks show that the above supremum is always achieved in 
R. So we can define t„ such that An = e"^^" |a^(t„)|. 

As in Step 2, we define t„ > to be the nearest local maximal point of the function 

Vp^it) to the point t„. We distinguish several cases according to the behavior of the 
sequence (in)n- We define the function by 

ai{t) = ^e-'*-ai{t + in). 

We can assume that, up to a subsequence, the sequence (a^)n converges on compacts 
to Ooo a nontrivial solution of 

o2 ^ (AT -2)2 N(N + 2)^ ,_2 

Of C^oo - Xj Ooo ^ Ooo H ^ (cosh t) Ooo = 

in the case where the sequence (t„ — in) is bounded, or to a nontrivial solution of 

2. ~ (iV--2)2^ 

C't ^00 Xj (X(x) ^ CLqq — 

in the case where the sequence (t„ — in) is unbounded. 

Moreover, Oqo is bounded by a constant times e^*. However, the choice of (5 G 
(—^^-2^,—^) implies that 6 < —~fj for all j = 0, ...,A'' and there are non nontrivial 
solutions of the above homogeneous problems which are bounded by e''* at +00. Hence, 
Coo = 0, a contradiction. This completes the proof of the result. ■ 
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Wc recall some well known result concerning harmonic extension of functions which 
are defined on S^~^. 

Lemma 8 Given ip G C^'°'{S^~^), we define to be the unique harmonic extension of 
(f in B (0,1), namely 



(105) 



AV^ = in 5(0,1) 

= (p on dB{0,l) 

Assume that ip is L'^{S^~^) orthogonal to eo, . . . , e^, then 

ll^¥'llcO(B(0,l)-{0}) ^ c||(^||cO(s]v-i) 
for som,e constant c > which does not depend on ^p. 

Using the fact that Kelvin's transform of an harmonic function V 

Wix) = \x?-^Vi^) 
\x\ 

is harmonic, the above result translates into the : 

Lemma 9 Given ip G C'^'°^{S^~^), we define W^p to be the unique harmonic extension of 
ip in IR^ — -B(0, 1) which decays at oo. Namely 

AW^ = in R^-B{0,1) 

(106) 

= (p on dB{0,l) 

and Wip tends to at oo. Assume that ip is L'^{S^~^) orthogonal to cq, • • • , ejv then 

\\^f\\c°_j^(R^ -3(0,1)} < c||(/?||cO(5iv-i) 
for some constant c > which does not depend on ip. 
Prom now on we assume that O is a bounded regular domain in IR^. 

8 Bubble tree solutions in general domains 

As before, we only prove the case when N > 6 since the proof of the result when N = 5 
follows the same lines with minor modifications. We recall 

2 

We define the space 

S := l^ipeC'^''^{S^~^) : J ^ ^ipejdio = 0, j = 0,...,N and \\ip\\c2,c < 
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8.1 Solution of the nonlinear problem in Qint,e' 

Given a m functions (f := {(pi, . . . , (pm) € and m points x := (xi, . . . , Xm) € we 
construct a positive solution of problem (68) in i^int,e whose boundary is, in some sense, 
parameterized by if. Namely we would like to solve 

^Uint,i + Xuint,i + U^riti = ^ B{xi,re) 

(107) 

Uint,i e Span {eo, . . . , Cat} on dB{xi, r^) 

For each i = l,...,m, we denote by the unique harmonic extension of ipi in 
B{xi,Re), namely 

AV;,^ = in B{xi,re) 

(108) 

V^^ = ipi on dB{xi,re) 

It follows from Lemma 8, together with a scaling argument, that 

Wfi\\e°[B{xi,re)-{xi}) < CrJ^ ll<^j|lcO(5JV-i)- (109) 

We keep the notations of the previous sections and, we look for a positive solution of 
problem (68) in B{xi,re) of the form 

Uint,i = UpXiii- - Xi) + V^. + Wi (110) 

where the function Up^x^^^ is the radial solution of problem (68) which has been obtained 
in Proposition 7 and where the functions Wi is small. 

As usual, we introduce the polar coordinates {t,9) G (— logr^j+oo) x S^~^ in each 
B{xi,re)- Given a function v, defined on B{xi,rc), we agree that the function v is the 
function defined on (— logr£,+oo) x S^~'^ which is determined by the relation 

2 

^(a;) = |a;| p-^ v{—log\x\,9). (HI) 
With these notations, we need to find a function Uint i 

and ■ ■ ■ 5 ^AT G K such that 

dfuint,i - ap dtUint,i - bp Uint,i + A^iv-iUint^i = - A e"^* Uint,i - U^nt,i (l^^) 

in [— logr£,+oo) x S'^~^ and 

2 N 

Uint,i{- log re, 9) = ri'^ ipi{9) + ^ bj ej 

j=o 

on S^-\ 

We will obtain a solution of this equation as a fixed point for some contraction mapping. 
We fix G (-(f + j^), -f ) such that {^ + 6- Tvq^)^ > 2 and we define 

£int,e ■.= i^weC'i{[- log r„+oo)xS^-^): ||ii;||co < Ke5+(^+'5-ivT2)ivb| (113) 
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where the parameter k> will be fixed later on. 
We write (112) as 

jrwi = -Xe-^\wi + %.)-Q^,{wi) (114) 
where the linear operator C is given by 

C:=dt + AgN-i - Qpdt-bp + p u^~x,^. 

and where collects the nonlinear terms 

Q^,{wi) := (up,A,€i + %i + - -Pul;^,^. Wi. 

We estimate ^ 

\\Xe-'^%,\\co < cA£5 ri^^^^^ (115) 

and 

\\Xe~'^^w\\c° < cXke^^'^^^^'^'^^ rl- (116) 

In view of the asymptotic expansion of UpjA.^i we have obtained in Proposition 7, it is easy 
to check that, for all w G Sint,e 

in (— logr£,+oo) x S^~^. Moreover, it follows from (109) that 

\V^i\ < cr-^ \\<Pi\\L-o e"^* < cr-^ e"^*, (117) 
in (— logr£,+oo) x S^"^. Hence, we conclude that 

Taylor's expansion yields 

(1 + *)'' - 1 - pt < ct^ 
near t = 0. This, together with the fact that b < —j^, implies that 

\\Q^,{w)\\co <c£^^^^^^^~^^^il + c^e^^-T^^), (119) 

5 

for some constant > depending on k. Wc have used the fact that + 6 — j^:2 ^ 1- 
Gathering the previous estimates, we conclude that 

\\-Xe-^\w + %,)-Q^,{w)Lo < c(1+c«£^)£^+(f^+*-^)^ (120) 

where > depends on k and the positive number 7 is independent of p. 
Given w G £int,e we use the result of Proposition 8 to solve 

It follows from Proprosition 8 and the above estimate that, given k, there exists £0 > 
(depending on k) such that the mapping 

'■ £int,e ^ £int,£ 
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defined by Ti{w) = v is well defined, provided e G (0,£:o)• 
Moreover, for all wi,W2 G £int,ei one can check that 

\\Ti{wi) - Ti{w2)\\cO < cX\\Wl-W2\\cO+c\\Q^^{Wl)-Q^.{w2)\\cO 
< c(A + £ 2 ) — tt;2 11(^0. 

S 

Consequently, for p sufficiently close to pi\r, the mapping Tj is a contraction from £int,e 
into itself and hence admits a unique fixed point in this set. This yields a solution Uint,i 
of (107). 

If we define the function Uint to be equal to Uint,i on B{xi,re)-, we have proven the : 

Proposition 9 Given x G Q,"^ and G S'^ , there exists a positive solution Uint of (68) 
in ^int,e satisfying boundary conditions 

Uint\dBix„r,) - ft ^ Span{ej : j = 0,...,N} 

for all 1 < i < m. Moreover, the sequence of solutions Uint blows up at each Xi as p tends 
to pn in such a way that 

m 

\VUint\'^ dx ^ ii 

i=l 

(3) 

in the sense of measures. Here C]y is the constant defined in Theorem 1. Finally, this 
solution can be expanded as 



1 



Clf^ (iV-2); VClf^ (2-N)!, {N-6)i 
6 2 + -I —e 2 X ^—^ ; P 2 



1 



+V^^ + 0{eirl) 
in B{xi,2re) - B{xi,rs/2). 

Since we have found the solution of (68) with the form (110), we have 

-Aw, = X{w, + V^,) + {up,x& +Wi + V^,r - ul^^^^^ (122) 
so that in B{xi, r^) — B{xi, re/2) 

\Awi\ < \wi + Vlp^\ + cu^~^^^.\wi + V^^\ < c\wi + 
Using the standard elliptic theory, we have 

Recall 

,2 ,2 ,A^ + 2 

(?3i+^ + *-ivT2)— >^ 



Thus, 

By the regularity theory, for all a G (0, 1 



-2 

reVtt;i||ioo(B(a;i,r,)-S(ii,3re/4)) < C£2 



rednWi\\ci,c(^S'^-'^) < c£2 (123) 
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8.2 Solutions of the nonlinear problem in Qext,e 

Given a m functions = (^i, . . . , 0^) G S'^, we now construct a family of positive solution 
of (68) in Qext,e which in some sense is parameterized by (f). 

Let X be a C°° cut-off function defined in R^, such that x|B(o,ro) = 1 ^'^d x = 
on — 5(0, 2ro) and x ^ 0- Denote by W^^ the unique harmonic extension of (pi in 
— B{xi,re) which decays at oo. We look for a solution of (68) in 0,ext,e of the form 



m 



^^ea;t = X! (^i^' ^(■'^») +^(" ~^») (^4>i + 
1=1 



ai-{- - Xi) 



+ w. 



ext 



(124) 



m 



where a := (ai, . . . , a^) G 
satisfy Wext\dne^t,e = 0- 

We use the maximum principle to reduce (68) to 



and the function Wext is assumed to be small and to 



-AWext = >^'Wext + Q + QA,<l>,aiWext) in ^ext,e 
Wext = on dQext,e 



where 



and where the function q is given by 



J2 l^i£~^G{; Xi) + x(- - Xi) (W^. + 



' (* Xi 



+ W 



q{z) = J2^x{z-Xi)(W^,{z) + 



i=l 



ai - {z - Xi) 



+ 2^Vx(^-x,)-VUy^,(z) + 



i=l 



Z-Xi\^ 

ai - {z- Xj) 



i=l 



+ \Y,(Kie^2G{z,Xi) + x{z-Xi)[W^.{z) + 



ai - {z - Xi) 



Given Aq and k > 0, we define 

g = {A G M™ : |A| < Ao} and = {a G {R^)'' : \a\ < rf }, 

Furthermore, given u e (2 — iV, 3 — AT) , we consider 

£ext,e = {w e Cl^{Qext,e) ■ WMlcS < r"^'" and w\dn,^t^^ = 0}, 
For all a G Ag, A G ^ and (pj G £, we estimate 

and given w G £ext,s, we obtain with little work 

\\Xw\\co , <c||Au;||co <c£^||u;||co 



(125) 



(126) 



(127) 



(128) 



(129) 



(130) 
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and 

\\QK,U^)\\ci_^<c{Kl + l + K^e'^)e'^r-^+^ (131) 
Finally, we estimate for all wi,W2 € £ext,£ 

2 

\\Q\,4,,a{:ilJl) - QA,4>,aiw2)\\cO_^ < CE^+i \\wi - W2\\c2 (^^^) 

The following result is standard 

Lemma 10 Assume that u ^ {2 — N,0) then for all f G C^-2{^ext,e), there exists w G 
C^i^ext,e) unique solution of 

Aw = f in ^ext,e 

(133) 

W = on d^lext,e ■ 

Furthermore, there holds 

\Mc2 ^ 4f\\co_^- 

Proof. The existence of w is straightforward and the estimate relies on the fact that 
X — |x — Xi\'^ can be used as a barrier in B{xi, ro) — B{xi, r^). ■ 

We define the map 

by <^ (j(u;) := V where v is the solution of 

Av = Xw + q + QK,ct>,a{w). 

Given k > 0, it follows from the estimates (129), (130) and (131) that the mapping 
TK,(t>,a is well defined and is a contraction, provided e is chosen small enough, say e G (0, eo). 
In particular, this mapping has a unique fixed point in £ext,£ which yields a solution of 
(125). Therefore, wc have proved the following : 

Proposition 10 Given x G $7"^, a G (M^)'" and (f) G S"^, there exists u^xt positive 
solution of equation (68) in ^ext,£, satisfying 

Uext = (Pi-\ N |-2^Aj£2 G{-,Xi) 

on dB(^Xi, Tg) for all \ ^ i ^ m and u^xt 

= on dU. Furthermore, the function Uext can 

be expanded as 

Uext = W^, + ' ^'^ '''' + ^ A,ei G{;Xi) + 0{erl) 

inB{xi,2re)-B{xi,re/2). 
Similarly, 

\\rednWext\\c^.c(s'^-i) <CS2r'^ (134) 

where n is the outside unit normal vector on the boundary of B{xi,re)- In the following 
consideration we will fix some a G (0, 1). 
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8.3 The cauchy data mapping 

We explain how the free parameters in Proposition 9 and Proposition 10 can be chosen 
so that the functions Uint,i and Ugxt can be glued together to obtain a positive solution of 
problem (68) in fl. 

We set ^ = (^1, . . . ,^m)- We want to choose the suitable parameters 

S := (x,A,(/j,^,a,0 

so that Uint.i and u^xt have the same Cauchy data on each dB{xi,rs)- Once this is done, 
the function defined hy u = Uint,i in B{xi,r£) and u = Ugxt in ^ext,e will be and solution 
of (68) away from the dB{xi,rs)- Elliptic regularity theory will then imply that it is a 
solution in J7. Moreover, it will follow from the construction itself that u has the desired 
behavior near each Xj and this will complete the proof of Theorem 1. 

Therefore, it remains to solve, for alH = 1, . . . , m, the system 



(135) 



9n'^int,i ~ ^n'^ext) 
on dB{xi,ri;). 

We denote by Ilj the L^(5"~^)-projection onto Spanjcj}, and 

N 

n(</,) :=</,-^n,(0) 

j=Q 

For alH = 1, . . . , m, the L^(S'"~^)-projection of (135) over the orthogonal complement of 
Span{eo, . . . , cat} yields the system of equations 

ipi = ^i + Fj,i(S), 

(136) 

Te d„ F^, = re + ,2 (S) , 

Next, we use the expansions of Lemma 1, Corollary 3 and Corollary 4 to obtain the 
L^(5'"'~^)-projection of (135) over Spanjeo} 

= e^Ai (r^-^ - H{xi, Xi)) + ^ AiG{xi,xi) + Fi,3(S), (137) 

Finally, the L^(5'"~^)-projection of (135) over Spanjei, . . . ,ejv} yields 

ai + s^lre^AiVzG{xi,xi)-r£AiVzH{xi,Xi)\ = Fi,5(- 



aiil- N) + e^ ireY^AiVzGixi,xi) -reAiVzH{xi,Xi)\ = Fi^^il 



(138) 
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Here for z = 1, . . . , m and I = 1, ... ,6 are continuous maps satisfying 

|F,,KH)|=0(e^r2). (139) 

We define "Dirichlet to Neumann map" for any 

5 : n(C2'"(5^-i)) n(c^'"(5^-i)) 

by 

where (resp. W^) is the harmonic extension in the ball .6(0, r^) (resp. in — i?(0, r^)) 
defined in Lemma 8 and Lemma 9. It is well known that S is an isomorphism [11] the 
norm of whose inverse does not depend on e. 

Hence, (136) (137) and (138) are equivalent to the following system 

= Gi,i(E;), 



^Kil^ci^Y^J ' (140) 

ai = rf-iGi,4(S), 

Va.,^^(x,A) = e-^r-iGi,5(H), 

VA,J^^(a;,A) = £-5Gi,6(S), 
where Gi^iiS) for all Z = 1, . . . , 6 and for alH = 1, . . . , m are continuous maps satisfying 



and 



(1) _ ^^-^(-^ ^(2) _ 



'^AT — 77^; — 2 ^ 

1 



N iA\ 2 



Moreover, elliptic regularity Theory shows that all Gi^ii^) are compact operators. 

Assume that (x'^,A'^) is a non degenerate critical point of JT^. In particular, this 
implies that dJ^^, evaluated at this point, is a local diffeomorphism from a neighborhood 
of (x'^,A'^) on a neighborhood of in ]R"*(^+i). Using this we can write formally the 
system (140) as 

S = $(S), 

We set & := — irr— log | — , | for alH = 1, . . . , m. We consider the set 

A = S((xO,AO),£i) X f X X 5(^0, £l) 

where ei is some fixed small positive number. It follows from the above analysis that 
<I> : ^ ^ ^ is a continuous compact map. According to Schauder fixed point theorem, $ 
has a fixed point in A. This completes the proof of Theorem 1. 

Remark If N = 4, we have the similar results. In this case, we take A = A^i^p^- 
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